Multifragment production in Au+Au at 35 MeV/u 
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\ Multifragment disintegration has been measured with a high efficiency detection sys- 

, tem for the reaction Au + Au at EjA = 35 MeV . From the event shape analysis and 

1^ , the comparison with the predictions of a many-body trajectories calculation the data, for 

0> , central collisions, are compatible with a fast emission from a unique fragment source. 



The disassembly of highly excited systems remains an open problem in the in- 
vestigation of intermediate energy Nucleus- Nucleus collisions []I|, One of the 
^ ■ challenging questions for head-on collisions is whether light particles and fragments 

emission is compatible with the fast emission from a unique thermalized source or 
^ I it can still be explained in the deep-inelastic framework. 

\ Several recent experimental studies of central collisions, performed with very 

heavy nuclei at different incident energies, give different indications on this 
point 0, ^, 1^, 1^. In 100 MeV /u Au + Au central collisions, dynamical and sta- 
tistical analyses suggest that the large multiplicities, observed for light particles 
and Intermediate Mass Fragments, are compatible with the prompt multifragmen- 
tation of a heavy, thermalized composite system with freeze-out density ~ | ^ | 
of the normal nuclear density (po = 0.15 fm~^), even if the fragment probability 
emission resulted strongly influenced by the radial flow |0, If. In the nearly sym- 
metric Pb + Au reaction at 29MeV/u [Q the charged products emission, studied 
for increasing neutron multiplicity, shows that the emission of Intermediate Mass 
Fragments becomes the largest component of the cross section at the expense of 
Projectile Like Fragments and Fission Fragments emission. On the other hand an 
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analysis, mainly based on the event shape, of the same reaction Pb + Au |^ at the 
same incident energy seems to reveal that, even selecting the most central collisions, 
the largest part of the total reaction cross section is due to strongly damped binary 
collisions. Other studies on light particle and fragment emission seem to confirm 
that, even selecting the most central collisions, the Incomplete Fusion cross section 
vanishes when reactions involving heavy projectile and targets at incident energies 
greater than ^ 35MeV/u are studied. 

In this Letter we report the results of the analysis of central collisions for the re- 
action Au + Au at 35 MeV/u, measured with a high efficiency detection system. We 
will show that the observed fragment emission is compatible with the fast emission 
from a unique equilibrated intermediate nuclear system. 

The experiment was performed at the National Superconducting Cyclotron Lab- 
oratory of the Michigan State University. Beams of Au ions at E/A = 35 MeV 
incident energy, accelerated by the K1200 cyclotron, were used to bombard Au foils 
of approximately 5 mg/cm^ areal density. Light charged particles and fragments 
with charge Z < 20 were detected at 23° < 9iab < 160° by 159 phoswich detector 
elements of the MSU Minihall [0. Reaction products with charge Z < 83 were 
detected at 3° < Oiah < 23° by the Multics array [|11[]. The charge identification 
thresholds were about 2, 3, 4 MeV/u in the Minihall ioi Z= 3, 10, 18, respectively 
and about 1.5 MeV/u in the Multics array independent of the fragment charge. The 
geometric acceptance of the combined apparata was greater than 87% of 47r. 

From the total charged particle multiplicity N^. the reduced impact parameter 

1) was determined, following Ref . [|l^ . Additional high statistics measurements were 
done using a shield, covering 9iab < 8° in order to minimize the radiation damage of 
the most forward detectors [|l^]. In this way more than 10^ events were collected for 
a centrality cut N,, > 24, b < 0.3, which selects 10% of the total measured reaction 
cross section. In this range of impact parameters the measured light particles {Z < 

2) multiplicity Mi has a mean value ~ 20 and the fragment {Z > 2) multiplicity 
Nf distribution has a gaussian shape with mean value 5.6 and standard deviation 
1.8. The mean value of M/ is very much higher than the one chosen in Ref. |^ 
{Ml > 10) to identify the central collisions. The requirement M/ > 10, corresponding 
to Nc > 15, would select an impact parameter range b < 0.6. Moreover the gaussian 
Nf distribution looks very different from the distribution shown in Fig. 2 of Ref. , 
where the two-fragments events represent the largest part of the measured cross 
section. 

To investigate the fragment emission patterns we first performed a shape anal- 
ysis, looking at the sphericity, coplanarity and fiow angle, variables sensitive to the 
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dynamics of the fragmentation process |T^. The emission of fragments from a unique 



source should be on the average isotropic in momentum space and the event shape 
should fluctuate around a sphere. Conversely in peripheral reactions the forward- 
backward emission of fragments from the spectator-like sources should lead to an 
event shape elongated along the beam axis, to a decrease of the sphericity value and 
to flow angles peaked in the forward direction. 

In this analysis only events satisfying the constraint that at least 70% of the in- 
coming momentum had been detected, were considered. For central and intermedi- 
ate impact parameters, where the particle and fragments detection is less influenced 
than in peripheral collisions by the energy thresholds and the angular acceptance, a 
further constraint on the total detected charge (70% of the total charge) was applied. 

The momentum tensor has been evaluated: 

(n) (n) 

T,=S^i-^ (.,J = 1,2,3) (1) 

where p'f^^pf^ are the i — th and j — th cartesian projections of the momentum p*^") 
of the n — th fragment in the centre of mass frame. The sum runs over the number 
of fragments {Z > 2) detected in each event. The diagonalization of the flow tensor 
gives three eigen- values and three eigen- vectors el. The event shape is an oriented 
ellipsoid with the principal axes parallel to the eigen-vectors. The sphericity and 
coplanarity variables are, respectively, defined as: 

5 = 1.5-(1-Ai), C=^-(A2-A3) (2) 

where Ai, A2, A3 are the ordered eigenvalues (Ai > A2 > A3), normalized to their 
sum. The flow angle 9fiow is the angle between the eigenvector el for the largest 
eigenvalue Ai and the beam axis: 



cos 



(Ofiow) = ei-k (3) 



Events with more than two detected fragments were used in the event shape analysis: 
two-body events, indeed, correspond mainly to peripheral reactions and do not give 
significant information, being two of the three eigenvalues zero for the momentum 
conservation. 

In Fig. 1 the S — C plot is shown, together with the cosine of the flow angle for 
three different gates on A^^- For < 15 which corresponds to 6 > 0.6, as expected 
for peripheral events, the memory of the entrance channel dominates: one or two of 
the three eigenvalues extracted from the momentum tensor are nearly 0, so that the 
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event is pencil or disk shaped and cos{6fiow) is peaked in the forward direction. With 
a rough constraint on the centrahty, i.e. requiring Nc larger than 15 (b < 0.6), the 
events do not show a well defined shape, though the centroid of the events is shifted, 
with respect to the more peripheral collisions, towards the corner which represents 
spherical events. This behaviour refiects in the cos{6fiow) distribution, which is 
less forward peaked than in the previous case. A more stringent constraint on the 
centrality {N^. > 24 , 6 < 0.3) leads to a drastic change of the S — C pattern: in these 
collisions the eigenvalues are very similar to each other, so that mainly events with 
shape close to a sphere are clearly present and cos{Ofiow) is randomly distributed, 
as expected in the case of fragments emitted from a unique source. Taking into 
account the modifications of the event shape with the requirement on the centrality 
and considering that the applied criterion on the total detected parallel momentum 
and the total charge does not eliminate the heavy residues, if they exist, we can 
deduce that the contribution from deep inelastic reactions is negligible at such small 
impact parameters |TB |. 

The next step of the shape analysis consisted in the comparison of the exper- 
imental mean values < S > and < C > and their standard deviations AS*, AC 



with the prediction of a many-body trajectory calculation |T^, which has as basic 
assumption the uniqueness of an emitting system with zero angular momentum. 

The experimental data considered in the following were selected with the central- 
ity cut Nc > 24 (6 < 0.3). In the simulation both charge and energy of the reaction 
products are selected by randomly sampling the experimental single-particle yield 
for this cut. The fragments are emitted from a spherical source of radius Rs and 
charge Zg. The individual emission time for each fragment is assumed to follow an 
exponential probability distribution, characterized by a decay constant r. The emis- 
sion is isotropic in the reference frame of the emitting system. A possible collective 
radial expansion can be accounted for by a further parameter Vcoiu which allows 
to increase the fragment velocity by a component v{r) = Vcoii^, which attains its 
maximum value at the surface i?, of the source. The simulated events were treated 



in the same way as the experimental data, after filtering |[T^ them with the geo- 
metrical acceptance, granularity, energy thresholds and finite energy and angular 
resolutions. In our case, due to the mass symmetry in the entrance channel, the 
fragment source is at rest in the centre of mass frame. This give the advantage that 
no hypothesis is needed on the source velocity (on the percentage of the momentum 
transfer) contrary to the case of asymmetric reactions |T^, when calculating the 
laboratory fragment velocities. 

The comparison between < 5* >, < C > is significant and represents a check of 
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the compatibility between the experimental observables and the decay of a unique 
source, provided that a set of input parameters {Zs, Rs, t, Vcoii) is found, reproducing 
the fragment momenta used in equation (1) to build the momentum tensor. To this 
aim we performed several calculations, varying in a wide range the input parameters. 
For each set the predicted fragment emission velocities as a function of the emission 
angle and the fragment reduced velocity correlation functions have been compared 
to the experimental data. The experimental N^., Nf and Zhound (total charge bound 
in form of fragments with Z > 2) spectra were continuously used to keep under 
control the reasonability of the predictions. 

We found that a source with charge ~ 86% of the total charge {Zs = 138), 
freeze-out density ~ Po/4 (radius Rg = 13/m), which emits the fragments with 
an average time between successive fragment emissions t ^ 85/m/c reproduces 
the previously mentioned observables, provided that an expansion radial velocity 
~ lAcm/ns is taken into account. As can be seen in Figs. 2 a), b) and c) the 
experimental distributions of the fragment emission velocities as a function of the 
emission angle 6cm are very well reproduced, irrespective of the selected fragment 
charge, even at forward and backward centre of mass angles, most affected by the 
experimental acceptance. 

Since the isotropy of the fragment emission, assumed by the calculation, implies 
that the emission velocity, for a given fragment charge, is constant over the whole 
range of 6cm, we investigated more in detail the main reasons of the rise of the 
experimental distribution at small centre of mass angles and of the dip at backward 
6cm- It is important indeed to understand whether these distortions can be ascribed 
to the boost of sources not at rest in the centre of mass or they are only due 
to experimental limitations. A simple kinematical calculation showed that, starting 
from a constant emission velocity with gaussian profile (in the calculation of Fig. 2 d) 
Scm/ns with standard deviation 1cm/ ns), the combined effects of the laboratory 
angular and velocity acceptance do not sharply cut the spectrum but they select 
at the most forward (backward) angles the highest (lowest) values of the velocity 
distribution. In particular the enhancement at 6cm < 50° is mainly due to the 
forward angular limitation and the dip at 6cm > 130° to the velocity thresholds. 
For 6cm ~ 50° ^ 130° the distribution is only slightly affected by the experimental 
limitations, so that the mean value in this angular range corresponds to the true 
emission velocity. For sake of comparison we report in Fig. 2 d) the experimental 
emission velocity for fragments with Z = 6, which have in the 6cm range 50° -r- 130° 
a mean emission velocity ~ 3cm/ns. 

The reproduction of the relative fragment momenta was then checked through 
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the comparison of the two-fragments correlation functions. The shape of the corre- 
lation function at small values of the reduced velocities is a measure of the spatial 
separation of the emitted fragments and is therefore sensitive to the input parame- 
ters of the calculation. 



The correlation functions ||19| , |2(| were calculated by: 



J- mix y'^red ) 

where Vred is the reduced velocity of fragments i and j {i ^ j) (charges Zi and Zj): 

I Vi - Vj I 

Vred = —P (5) 



Y{vred) and Y.mix{vred) are the coincidence and mixed yields for fragment pairs of 
reduced velocity Vred- The mixed yield was constructed by means of the mixing 
event technique, C is a normalization factor fixed by the requirement to have the 



same number of true and mixed pairs [19 



We analyzed separately fragments detected in the Multics array and in the Mini- 
hall, since the solid angle covered by the apparata is very different and an average on 
the whole solid angle would lead to a loss of information. From Figs. 3 a), b) (frag- 
ments detected in Multics) and Figs. 3 c), d) (fragments detected in the Minihalt) 
one can see that the many-body trajectories code well reproduces the experimental 
correlation functions, irrespective of the selected fragment charge, assuming an aver- 
age time between successive fragment emissions r = 85 fm/c and a collective radial 
expansion Vcoii ~ lAcm/ns. Varying r by some tens of fm/c the experimental cor- 
relation functions are not so well reproduced: the decreased or increased distance 
among the emitted fragments introduces additional correlations or anticorrelations. 



respectively, not present in the data [21]. It has to be noted that the experimentally 
observed Coulomb hole at small values of the reduced velocity leads to the same 
choice of the parameter r both in the case of fragments emitted at small relative 
angles (detected in Multics) and in the case of fragments emitted at large relative an- 
gles (detected in the Minihall). In addition in the first case the selection of r can be 
performed by checking even the reproduction of the enhancement of the correlation 
function at Vred ~ 20. This bump, due to the mutual Coulomb repulsion between 
the closely emitted partners (see Fig.s 3 a) and b)), is sensitive to the increase of 
the proximity of the fragments induced by the decrease of r. 

The small value found for the radial collective velocity, although consistent with 
the extrapolation to 35MeV^/ u of data at higher energies [^ , should be thoroughly 
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investigated before making any conclusion. It could be ascribed either to the as- 
sumption of non overlapping fragments or to the schematic treatment of the Coulomb 
component: a change of the emission geometry, which reflects in an increase of the 
fragment Coulomb energy, could compensate the need of a radial expansion. Fur- 
thermore in a recent theoretical work on the effect of the angular momentum on 
the statistical fragmentation, it was found that a rotation mechanism could explain 
some features previously ascribed to a collective flow of the nuclear matter. The 
comparison to predictions of statistical models which take into account either the 
radial flow or the angular momentum for systems with mass of some hundreds of 
nucleons should be performed, but this is beyond the aim of this Letter. 

Since the observed fragment emission patterns are well reproduced by the many- 
body trajectory calculation, the comparison between experimental and predicted 

< 5* >, < C > and their standard deviations AS*, AC as a function of the fragment 
multiplicity becomes significant to draw a conclusion on the compatibility of our data 
with the uniqueness of the decaying system. From Fig. 4 the very good agreement 
between data and predictions is evident, not only for < S' >, < C >, but even for 
the standard deviations AS, AC, irrespective of the selected fragment multiplicity. 
Observing the unfiltered predictions (dashed lines in Fig. 4), we can deduce that the 
effects of the experimental inefficiencies very slightly decrease < 5* > and increase 

< C >. Moreover events with Nf = 3, not present in the unfiltered predictions, 
are only few percent both in the experimental data and in the filtered predictions. 
We would like to stress that this agreement is not a trivial consequence of the 
reproduction of the fragment multiplicity. Indeed even in the case of peripheral 
collisions (Fig. 1 a)) events with fragment multiplicity up to 7 were measured, but 
in this case the forward-backward emission flattens the fragment momenta into a 
plane and < S > and < C > reveal pencil/disk shaped events. 

In conclusion, for the central collisions of the reaction Au+Au ai E/A = 35 MeV, 
the good agreement among the measured observables and the predictions of a many- 
body trajectories code confirms both the assumptions on the uniqueness of the 
fragment source and on the isotropy of the fragment emission. 

This work has been supported in part by funds of the Italian Ministry of Uni- 
versity and Scientific Research. The technical assistance of R. Bassini, C. Boiano, 
S. Brambilla, G. Busacchi, A. Cortesi, M. Malatesta and R. Scardaoni during the 
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Figure captions 

Fig. 1: Experimental Sphericity- Coplanarity linear contour plot and cos{9fiow)- Upper 
panels: Nc < 15 (6 > 0.6), intermediate panels: Nc > 15 (6 < 0.6), lower panels: Nc > 24 
(6 < 0.3). 

Fig. 2: a), b), c): Experimental emission velocities (points) for fragments with charge 
4, 8 and 10, compared with the many-body trajectory calculations (line). The statistical 
experimental error is smaller than the size of the points. 

d) Experimental emission velocities for fragments with charge 6 (open points) compared 
with a constant emission velocity (3cm/ns with standard deviation Icm/ns) (line) fil- 
tered with the constraints 8° < Oiab < 160", Vthreshold = l-5cm/ns {9iab < , Vthreshold = 
2.5cm/ns {Oiab > 23°). 

Fig. 3: Two fragment correlation functions for: a) 3 < Zjmf < 30 and 8° < Oiab < 23°, 
b) 3 < ZiMF < 20 and 8° < Oiab < 23°, c) 3 < Zjmf < 20 and 23° < Oiab < 160°, 
d) 3 < ZjMF < 10 and 23° < Oiab ^ 160°. Open points show experimental data. 
The solid, dashed and dotted lines are the the many-body trajectory predictions for 
r = 85, 50, 150/m/c respectively. 

Fig. 4: Mean sphericity < S > and coplanarity < C > and their standard deviations 
AS, AC as a function of Nf. Experimental < S >, < C >,AS,AC are reported as 
full symbols and vertical bars, respectively. < S >,< C > calculated from the filtered 
predictions are reported as open symbols and < S > ±AS, < C > ±AC as full lines. 
< S > ±AS, < C > ±AC from the unfiltered predictions are shown by dashed lines. 
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